Intrinsic Spin Hall Edges 
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The prediction of intrinsic spin Hall currents by Murakami et al. and Sinova et al. raised many 
questions about methods of detection and the effect of disorder. We focus on a contact between a 
Rashba type spin orbit coupled region with a normal two-dimensional electron gas and show that 
the spin Hall currents, though vanishing in the bulk of the sample, can be recovered from the edges. 
We also show that the current induced spin accumulation in the spin orbit coupled system diffuses 
into the normal region and contributes to the spin current in the leads. 

PACS numbers: 72.25.Dc, 72.25.Mk, 72.20.Dp 
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Transport and manipulation of spins in semiconduc- 
tor structures has become a mainstream in condensed 
matter physics Q. In principle, spins can be injected 
into semiconductors by ferromagnets via electric con- 
tacts. However, finding suitable material combinations 
that do not suffer from the conductance mismatch 0, 
turned out to be difficult. Furthermore, introducing fer- 
romagnetic materials into the semiconductor microfabri- 
cation process is undesirable from a technological point 
of view. The prospect to generate spin accumulations 
in semiconductors without ferromagnets or applied mag- 
netic fields simply by driving a current through a mate- 
rial with intrinsic spin-orbit (SO) interaction and broken 
inversion-symmetry 0- Q is therefore very attractive. 
A related effect that attracted a lot of attention is the 
spin Hall effect (SHE), i.e. the spin current(SC) that has 
been predicted to flow normal to an applied electric cur- 
rent in the absence of an applied magnetic field. When 
caused by impurities with spin-orbit scattering 0, Q this 
effect is called "extrinsic" . A spin Hall current (SHC) can 
also be generated by the spin-orbit interaction of the lat- 
tice potential as has recently been predicted for p-doped 
III-V semiconductors || and the two-dimensional elec- 
tron gas with a Rashba-type SO interaction (R2DEG) [9|- 
Whether the experimental observations of the spin Hall 
effect by optical methods |l(J have intrinsic or extrin- 
sic origin is still a matter of debate. In spite of initial 
controversies, analytic theorie s Fill llil Hi! 0] as 
well as numerical simulations |l6l Il7j consistently pre- 
dict that the SHE should vanish in the disordered (bulk) 
R2DEG 0. Some doubts remain whether the SC, be- 
ing a non-conserved quantity in SO coupled systems, is 
observable at all [lflj . 

In this Letter, we focus on the spin currents near nor- 
mal contacts. First, an elementary and general proof 
is given that the spin Hall effect due to the lattice SO 
coupling (viz. intrinsic SHE) must vanish in diffuse bulk 
systems with an arbitrarily strong SO interaction that is 
linear in the electron wave vector. Nevertheless, using 
an extension of this argument to finite system sizes, we 
show that near the edges a spin Hall current can persist. 
Next, by solving the kinetic equations for a model sys- 



tem of a R2DEG in contact with a normal metal system 
without SO interaction, we calculate indeed a finite SHC. 
This SC is generated in a skin depth determined by the 
Dyakonov-Perel Q spin- flip diffusion length (L s ) and the 
polarization is not normal to the 2DEG, having a compo- 
nent due to the diffusion current from the SO-generated 
spin accumulation (SA). The magnitude of the SC gen- 
erated at the edges depends on whether the system is 
clean (impurity broadening less than the SO splitting) 
or dirty (opposite limit). However, in contrast to the 
bulk SC, the edge SC does not vanish when the system is 
not ballistic (L s smaller than the system size). The SC 
is calculated in the normal metal contact and therefore 
certainly a transport current |19| • Related work on inter- 
face and boundary effects focused so far on mesoscopic 
systems via numerical simulations |2oj | and the SA near 
hard wall boundaries |2lj . 

We proceed to derive a transport equation valid in the 
Boltzmann limit that is capable of handling the full spin 
dynamics. In 2 x 2 spin space, the Hamiltonian is 



H = p 2 /2m + V(x) +H R - eE(t) ■ x, 



(1) 



where x and p are the (two-dimensional) position and 
momentum operators, respectively. Here the unit vector 
z is normal to the 2DEG, H H = (a/H)p ■ (a x z) is the 
Rashba Hamiltonian with Pauli matrices <x and a param- 
eterizes the strength of the SO interaction |2^[. E is the 

electric field, and V(x) = Yli=i <K X — -^i) i s tne impu- 
rity potential, modelled by N impurity centers located 
at points {X{\. Although it is possible to consider ac 
fields, we focus here on dc fields in the x— direction and 
assume that the electric field is turned on adiabatically 
in the remote past at which the system was in thermal 
equilibrium, i.e. we assume E = lim s ^o Ea exp(st/h). 
To leading order in the impurity potential the diagonal 
elements of the density matrix in reciprocal space satisfy 
the following equation [j^ : 



■isf(k) + [f(k), H?] = (fkk>V k , k - V kk ,f k , k 



eE-[x,f }. 



(2) 
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Here / = (/° + /° )/2 + <r e (/° - /° )/2 is the equilibrium 
density matrix with f±(k) — T(h 2 k 2 /2m± akj, where 
T{E) is the Fermi function and ag = k • (er x z)/k. The 
off-diagonal elements of the density matrix read 

fkh> = ^ [°° dEG%(E)(f(k) - f(k'))G%(E)V kk ,. 

(3) 

Here G* {A) {E) = {E-H -H^+{-)is/2)- 1 are retarded 
(advanced) matrix Green functions. Substituting Eq. @ 
into Eq. J3J) and averaging over V k k' (in the Boltzmann 
limit averaging is equivalent to replacing |Vfcfc'| 2 with its 
average value N\<f> kk i \ 2 /A 2 , where N is the number of 
impurities and A is the area 23] ) gives our basic equation, 
valid for weak V and low enough impurity densities to 
ignore weak localization effects, but to all orders in a. 

The mechanism behind the intrinsic SHE is the spin 
precession of quasiparticles while being accelerated by 
the electric field |9§. However, impurity scattering pro- 
vides a brake that in the steady state cancels the accel- 
eration on average. Therefore the SHE should vanish in 
an infinite, homogeneously disordered system. This idea 
can be formally expressed by considering the acceleration 
operator Xj = eEi/m~ "S/iV/m — 2a 2 h~^t 3 jiU3'Pj, where 
eijk is the antisymmetric tensor and the Einstein sum- 
mation convention is implied. We notice that the last 
term is proportional to the j'th component of the SC op- 
erator polarized in the z— direction, JJ = {vj,cr z }. The 
expectation value is defined by (O) = Tr/O, where the 
trace is over wave vector and spin space, and 7TT denotes 
averaging with respect to impurity configurations. In a 
steady state the average acceleration (x^) must vanish, 
leading to the equality 



2a 2 m 2 h" 



e 3ji 



<V,V>. 



(4) 



We show that the right hand side of this equality also 
vanishes by evaluating the expectation value of the de- 
celeration due to impurity scattering: 



(ViV) = - K)v kk ,tvf k>k 

kk' 

= -i$^tr (k i [f(k),H k H ]+ik i eE j V kj f)=eE i , 

k 

where tr is the trace over spin components and Eqs. J2J- 
|3J) have been used in the second step. Substituting the 
expression above into Eq. Q we see that all components 
of the SC polarized in the z— direction vanish with the 
average acceleration [24|. This result holds for infinite 
systems regardless of the range of the impurity potential 
or whether the system is clean (akpr/H ^> 1) or dirty 
(akpr/h -C 1), where r is the momentum lifetime. Thus 
generalizing previous results 0, 0, 0, Q, 0] . How- 
ever, as we shall show below, for semi-infinite and finite 
systems, SCs persist near the edges, but the size of these 
currents depend on whether the system is clean or dirty. 



This line of argument allows one to check related 
Hamiltonians. In the presence of fc-linear Dresselhaus 
and Rashba terms, the result remains unchanged besides 



the substitution a 2 



-0 1 , where j3 is the Dresselhaus 



spin orbit coupling constant. Thus the SHC still van- 
ishes (with the possible exception of the degeneracy point 
a = (3 25]). When the SO coupling contains cubic terms 
like a(k) — a + «ifc 2 , it is easy to show that the SHC 
is proportional to u\ [l^. Another possible situation is 
the presence of a Zeeman field: in this case the operator 
equation is modified to give ae 3 ji(Jj) = (a^B^ — {cr^)Bi, 
relating the SHC to the SA. If a varies in space, the SHC 
is found to be proportional to the spatial derivatives of a 
and /. If a is constant, but / varies, e.g. due to bound- 
aries or interfaces the SHC is proportional to the 
gradients of the density matrix: 



4a 2 rn 2 

h 3 



£3i<<J/) = (H -+a(zx ( r) h V l f(k,x)}). (5) 
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This equation shows that although the bulk SH cur- 
rent vanishes, there is no a priori reason for SH cur- 
rents near the edges of the R2DEG to vanish. Next, 
we shall show indeed the SH currents do not vanish near 
the edges. We therefore return to the quantum trans- 
port equation, allowing for spatially varying density ma- 
trices but assuming short-range s-wave scatterers with 
|Vfcfc'| 2 = NX 2 1 A and a/k F <C 1. We solve the transport 
equation by expressing / in terms of a gradient expansion 
of p(E) = ^h 2 /2nm)j: k (G*iE)f(k) - f(k)G£(E)). 
In the case of s-wave scatterers and to leading order in 
ma/h 2 kF, this generates the same diffusion equation as 
Ref. 0| and Burkov et al. in Ref. 0. In terms of com- 
ponents of the density matrix: p = n + s • <x + 03S3 : 



DV 2 n - 4X S _ C (V x s) z = 



(6) 



2s 3 



DV's 3 - 2K p (V ■ s) = — ^ (7) 



£>V 2 s + 2K p Vs 3 - K s _ c (z x V)n = — 



(8) 



Here D = v 2 f t/2, t s = r(l + 4£ 2 )/2£ 2 , K s _ c = < 2 /(l + 
4£ 2 ), K p = hk F £/m(l + 4£ 2 ) 2 and f = ak F r/h. The SC 
is given, in the diffuse limit by 



h = 



l + 4<£< 



(s l3 (t 



- $ijS3j —DVjSi. 

(9) 

Electric field dependence can be reintroduced by the sub- 
stitution V — ► V + eEd E - 

We now focus on a four terminal structure as depicted 
in Fig.^ This structure consists of two massive reservoirs 
biased to produce a charge current in the a;— direction. 
Between the reservoirs there is a R2DEG hybrid struc- 
ture, with a(x) = «o for L > y > and a(x) = for 
y < and y > L, and additional differences are dis- 
regarded. The normal 2DEGs are coupled to massive 
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FIG. 1: The schematic setup for SC generation in an R2DEG 
that is current biased by reservoirs Rl and R2 and equipped 
with 2DEG Hall contacts with vanishing spin orbit coupling 
to reservoirs R3 and R4. The latter can be magnetic or non- 
magnetic and voltage biased such that no charge current flows 
through the 2DEGs. The dashed arrows indicate the local 
SHC densities that are concentrated near the interfaces. The 
dominant part of the spin current flowing into R3 (or R4) 
is generated within a skin depth L a near the corresponding 
interface and contains a diffusion term from the SA S y in 
the bulk of the R2DEG. The contribution of the SHC density 
near Rl and R2 to the net SC [T3l flowing into R3 (or R4) is 
exponentially small. 

reservoirs R3 and R4, that are biased such that the charge 
current is zero, but a SC can still be collected. To the or- 
der (in a/kp) that we are considering, n does not depend 
on y. Alternatively, one can assume that the transverse 
size of the leads to R3 and R4 are much smaller than the 
distance between Rl and R2, in which case one can also 
neglect the y dependence of n. We consider the distribu- 
tions at a safe distance from the reservoirs Rl and R2. 
In this region s v and s z depend only on y and s x = and 
the diffusion equation becomes: 

where r\ = (1 + 4£ 2 )~ 3 / 2 , y = y/L s and L s = \JDt s . 

In order to derive the matching condition for the spin 
and charge distribution functions at the contacts (i.e. 
interface between the R2DEG and 2DEG), short-range 
fluctuations of boundaries and interfaces that can lead 
to additional spin relaxation [3(| are disregarded. We 
consider an arbitrary solution \ °f the Schrodinger equa- 
tion set by the Hamiltonian Eq. QJ. We label the 
solutions in the R2DEG and 2DEG regions \R ancl 
Xn, respectively. At the interface Xfllo = Xn\o and 
n- (iV + ao(z x <j))xr\o — * n ' ^Xn\o, where n is the unit 
vector normal to the interface. Multiplying from the left 
with x\ o &i and evaluating the imaginary part we obtain: 

n ' (ixlr^XN - i(Vxjv) t cr J XA') 

= n ■ {ix R ^XR - i^XR^ViXR 

+ aoX R Wi,(z x &)}xr) ■ 



^ S/xK^. r ± 1 

« dx ' 
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FIG. 2: SCs and accumulations in the R2DEG (y < 0) as a 
function of distance from the interface to the 2DEG. The full 
line is the local value of the SHC density for akFr/h — 0.1. 
At the boundary the SHC density recovers its maximum value 
jb = eEK s - c /2na (it approaches the universal value eE/8n 
in the clean limit £ 3> 1) • The dashed line is the 2-component 
of the diffusion current density. Inset: Corresponding local 
SA in the same system normalized to the magnitude of the 
bulk SA eEaTm/2TT. The solid and dashed lines represent 
the y and z components, respectively. 

We identify the right (left) hand side of this equation as 
the SC density in the Rashba (normal) 2DEG. In terms 
of spin density matrices we have ti(f R {ai, n • j(0)}) = 
tr(/ JV {tr i , n-j(O)}), where j(x) = {v, <5(x— x)} is the local 
current density operator. We therefore have to match the 
normal components of the SC density given by Eq. © at 
the interface Since the operator {v^, a-j} can have a 
nonzero expectation value in the equilibrium state it has 
been questioned whether it governs transport of spins 
in the presence of SO interaction [J^. We notice that 
the negative energy solutions (relative to the band cross- 
ing) of the Hamiltonian Eq. Q without the electric field 
term, are localized to the R2DEG region if surrounded 
by a region with a = 0. In a normal 2DEG surround- 
ing the R2DEG, we can therefore show that equilibrium 
SCs exposed in Ref. [lflj l do not transmit into the normal 
region. Moreover, the expectation value of the SC den- 
sity operator vanishes for these localized solutions and it 
is precisely the absence of contributions from these solu- 
tions that shifts the equilibrium value of the SC to zero. 

Returning to the setup in Fig. ^ we assume that the 
reservoirs R3 and R4 are sufficiently large such that all 
components of the SA at their respective interfaces with 
the ordinary 2DEG leads vanish. Shrinking the widths 
of the 2DEGs to zero, we obtain the effective boundary 
conditions Si — at the R2DEG|R3 interface. The fi- 
nite Ohmic resistance of a finite 2DEG region between 
the R2DEG and the reservoir can easily be reintroduced 
if necessary and would lead to somewhat smaller spin 
conductances. We then can solve the diffusion equation 
above and obtain the spin current using Eq. ©. 

Analytical formulae turn out to be too lengthy to re- 
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produce here. Our results are therefore summarized in 
Fig. |2I The SA is suppressed at the interface, reflecting 
the massive-reservoir boundary condition. The gradient 
of the two components s y and s z leads to two SC com- 
ponents. The SC polarized in the y-direction represents 
the out-diffusion of the bulk s y SA. This is not a Hall 
current, since it flows into the side contacts with oppo- 
site directions (Fig. EJ) with polarization that is inverted 
with the bias current direction. The resulting spin con- 
ductivity at the interface is a% y = 0.87e£ 2 /27r in the dirty 
limit (£ <C 1). For larger values of £, a v xy increases above 
this quadratic behavior, but in the clean limit (£ ^> 1) 
this increase is cut off by the resistance of the normal 
region. When spins diffuse from a finite distance into the 
2DEG, they precess in the SO-generated magnetic field. 
Consequently there is a diffusion current polarized along 
the z-direction, for which we find in the dirty limit a con- 
ductivity a z xy — 0.83e£ 2 /27r. The conductivity a*, con- 
trary to (jVy, decreases below this quadratic behavior for 
larger values of £ and vanishes in the clean limit. Nikolic 
et al. j2oj | recently observed SCs with z and y polariza- 
tion in numerical simulations. In addition, we also find 
a SHC exponentially localized to the edges that decays 
in the bulk on the length scale L s and reaches its maxi- 
mum value eEK s - c /2-Ka at the interface to the reservoir. 
This is due to the fact that the first term in Eq. be- 
ing proportional to Sj,(thus zero at the interface) is no 
longer screening the second term (proportional to Vn) 
and reflects the physical process that the SHC density 
generated near the interface can escape into the reservoir 
before it decays due to spin relaxation. The resulting dc 
spin Hall conductivity is given by e£ 2 /2ir in the dirty and 
e/8n in the clean limit. Our result differs from that of 
Ref. 01 w ho did not take into account the edge currents 
and obtained similar values only for the ac response at 
carefully tuned frequencies. 



In conclusion, we find that in a Hall geometry two 
different spin currents can be extracted by the Hall con- 
tacts from the current-biased disordered R2DEG. In ad- 
dition to the SHC, the current-induced SA drives a spin- 
diffusion current. The SO generated spin accumulation is 
therefore not confined to the region where it is generated, 
but can be extracted and, at least in principle, used as 
a source of spins for spintronics applications. Both dif- 
fusion and SHCs are generated within a strip that scales 
like the Dyakonov-Perel spin diffusion length. 
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